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Abstract: Understanding the effects of degradation on belowground biomass (BGB) is essential for 
assessment of carbon budget of the alpine meadow ecosystem on the Tibetan Plateau, China. This 
ecosystem has been undergoing serious degradation owing to climate change and anthropogenic activities. 
This study examined the response of the vertical distribution of plant BGB to degradation and explored 
the underlying mechanisms in an alpine meadow on the Tibetan Plateau. A field survey was conducted in 
an alpine meadow with seven sequential degrees of degradation in the Zoige Plateau on the Tibetan 
Plateau during the peak growing season of 2018. We measured aboveground biomass (AGB), BGB, soil 
water content (SWC), soil bulk density (SBD), soil compaction (SCOM), soil organic carbon (SOC), soil 
total nitrogen (STN), soil total phosphorus (STP), soil available nitrogen (SAN), and soil available 
phosphorus (STP) in the 0-30 cm soil layers. Our results show that degradation dramatically decreased the 
BGB in the 0-10 cm soil layer (ВСВо-10) but slightly increased the subsoil BGB. The main reason may be 
that the physical-chemical properties of surface soil were more sensitive to degradation than those of 
subsoil, as indicated by the remarked positive associations of the trade-off value of BGBo-10 with SWC, 
SCOM, SOC, STN, SAN, and STP, as well as the negative correlation between the trade-off value of 
ВСВо ло and SBD in the soil layer of 0-10 cm. In addition, an increase in the proportion of forbs with 
increasing degradation degree directly affected the BGB vertical distribution. The findings suggest that the 
decrease in the trade-off value of ВСВо-10 in response to degradation might be an adaptive strategy for 
the degradation-induced drought and infertile soil conditions. This study can provide theoretical support 
for assessing the effects of degradation on the carbon budget and sustainable development in the alpine 
meadow ecosystem on the Tibetan Plateau as well as other similar ecosystems in the world. 
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1 Introduction 


Grasslands occupy over 25% of the Earth's terrestrial surface and store approximately 30% of the 
global terrestrial carbon (C) pool (Hewins et al., 2015). Soil C sequestration in grasslands occurs 
at an estimated rate of 0.5 Pg C/a globally and accounts for approximately a quarter of the world's 
potentially sequestered soils (Qiu et al., 2013). The Tibetan Plateau, approximately 35% of which 
is occupied by alpine meadows (Zheng et al., 2000), is equal to about 44% of China's total 
grasslands and approximately 6% of grasslands worldwide (Piao et al., 2012). Alpine grasslands 
are vital to local livestock husbandry development and ecological functions, including soil and 
water conservation, biodiversity maintenance, and C sequestration (Li et al., 2013; Qin et al., 
2018). However, this region's grasslands are very fragile owing to the unique climate and harsh 
environmental conditions (Sun and Wang, 2016; Sun et al., 2019a). These grasslands have 
undergone serious degradation under the influence of climate change and anthropogenic activities 
(Harris, 2010; Sun et al., 2014; Dong et al., 2020), greatly threatening local economic 
development and ecological security (Huo et al., 2013). 

As a main source of the soil C pool and plant C stock, belowground biomass (BGB) plays an 
important role in the C cycles of terrestrial ecosystems (Mokany et al., 2010; Sun et al., 2013), 
especially in grassland ecosystems where over 25% of the plant C is stored in roots (Li et al., 
2008; Yang et al., 2009). On the Tibetan Plateau, where BGB constitutes the core portion of the 
grassland's total primary productivity (Li et al., 2011), the root-to-shoot ratio can be as high as 5.8 
(Yang et al., 2009). The alpine grassland ecosystem of the Tibetan Plateau is currently considered 
a C sink (Yang et al., 2008), largely because the roots account for a notably great proportion of 
alpine plants (Wu et al., 2011). Therefore, it is important to precisely estimate the root biomass 
and turnover rate to obtain a more reliable scientific assessment of the C budget in this terrestrial 
ecosystem. However, because of methodological difficulties and high cost, the distribution and 
transfer of a plant's ВСВ in response to grassland degradation is not well understood (Sun et al., 
2018a). 

Roots serve as both the primary pathways for a plant's nutrient and water uptake (Wu et al., 
2011) and its anchorage to the basement soil. As they consume resources, roots are involved in 
underground competition with other plants (Hortal et al., 2017), complex mutualistic relationships 
with beneficial bacteria and fungi (Duran et al., 2018), and competitive correlations with 
detrimental microbes (Rey and Schornack, 2013). The root system's overall structure and ability to 
adapt to habitat conditions is called plasticity and is a crucial determinant of the overall plant 
adaptability (Gruber et al., 2013; Rogers and Benfey, 2015). For example, the roots of perennial 
plants in particular legumes can affect the physical properties of forcefully compacted soils (Glab 
and Szewczyk, 2014). Because they enhance the water retention properties, these types of roots 
can be used to regenerate soils that have undergone physical degradation. However, some species 
have different effects on the soil environment, owing mainly to their unique root morphologies 
(Silva et al., 2014). In addition, root systems affect plant growth through nutrient foraging, gaseous 
exchange, and absorption of beneficial microorganisms (Pieterse et al., 2014; Verbon and 
Liberman, 2016). 

The main properties of the root system are as follows: length, biomass, volume, anatomy, and 
surface area (Lynch and Brown, 2012). These characteristics facilitate the ability of plants to 
survive under various biological and abiotic stresses, including salinity, extreme temperatures, 
waterlogging, drought, and nutrient deficiency (Dukes et al., 2005; Rodrigues et al., 2015). Thus, 
roots are considered to be an important factor determining the competitive ability of individual 
plants under environmental stress, which affects the responses of the plant community structure to 
environmental changes (Haichar et al., 2008). The characteristics and morphology of grass roots 
are closely linked to environmental conditions, and these roots tend to have highly dynamic 
lifespans ranging from a few days to weeks (Tingey et al., 2010). In addition, previous researches 
have revealed that the distribution of plant root biomass can be affected by various environmental 
factors, including soil physical-chemical and biological properties, climatic conditions, and plant 
characteristics (Rodriguez et al., 2007; Macinnis-Ng et al., 2010; Li et al., 2011). For instance, 
Zhang et al. (2019) found that plants allocate more root biomass to deeper soils for maintaining 
plant water acquisition and growth in response to drought stress. As grassland degradation 
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induces dynamics in the soil water, nutrients, microenvironment, and plant properties, the BGB 
distribution patterns might also fluctuate during degradation. There is little available information 
on the variation of the direction and magnitude of BGB distribution patterns in response to 
degradation of alpine meadows. We hypothesized that BGB in topsoil but not subsoil may be 
more vulnerable to degradation because surface soil loses water and nutrients relatively quickly 
compared with subsoil. 

In this study, we examined variations in BGB distribution patterns along degradation gradients 
in an alpine meadow on the Tibetan Plateau of China. The aims of the study were: (i) to verify our 
hypothesis on the BGB dynamics at different soil depths in response to grassland degradation; 
and (ii) to explore the underlying mechanisms driving the changes in BGB vertical distribution in 
a degraded alpine meadow. Our results may offer insights into grassland degradation and, more 
specifically, the effects of degradation on the C budget of the alpine meadow ecosystem. 


2 Materials and methods 


2.1 Study area 


The study was conducted in an alpine meadow at the Zoige Plateau (32°20'—-34°00'N, 101°30’— 
103930Е), which is located in the northeastern part of the Tibetan Plateau of China (Fig. 1). The 
Zoige Plateau has an average altitude of 3500 m and covers an area of 6180 km? (Liu et al., 
2020a). The climate is characterized as cold and humid. The mean temperature during the year, in 
the hottest month (July), and in the coldest month (January) is 1.5°C, 11.1°C, and —9.7°C, 
respectively (Liu et al., 2020b). The mean annual precipitation is approximately 700-800 mm; 
most precipitation occurs from May to August (Liu et al., 2020a). The plant community at our 
study site is dominated by Kobresia tibetica, Stipa capillata, and Carex lasiocarpa (Liu et al., 
2020b). The soils are classified mainly as Cambisols using the World Reference Base for Soil 
Resources (WRB) soil classification system and as Inceptisols using the United States 
Department of Agriculture (USDA) soil taxonomy (Liu et al., 2020a). 
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Fig. 1 Location of sampling site at the Zoige Plateau on the Tibetan Plateau of China (a), photograph showing 
the landscape of sampling site (b), and schematic illustration of experimental design (c) 


ZHANG Zhenchao et al.: Degradation leads to dramatic decrease in topsoil but not subsoil root biomass... 


2.2 Experimental design and soil sampling 


We sampled nine degraded grasslands and chose seven sequential degradation gradients according 
to the vegetation coverage and plant community structure, which is a common method of 
classifying degradation in this type of study (Ma et al., 2002; Zhang et al., 2019; Liu et al., 
2020a). All these grasslands were freely grazed; the livestock was mainly yaks and horses, which 
contributed greatly to the degradation. Table 1 displays detailed information on the seven 
degradation gradients, which were adjacent and thus had similar topography and environmental 
status; the heterogeneity within a degradation gradient was relatively low. The sampling survey 
was carried out in August (the peak growing season) 2018. We selected three representative 
sampling plots which were approximately 50 mx20 m as replicates on each degradation gradient. 
Next, three quadrats (1 mx1 m) within each sampling plot were randomly selected to obtain plant 
and soil samples for measuring their characteristics. The plant species were recorded, and the soil 
compaction (SCOM) was determined via a portable time-domain reflectometry equipment (TDR 
100, Spectrum Technologies Inc., Chicago, IL, USA). 


Table 1 Information on the seven degradation gradients 


Degradation 5 с Е : o 
gadiem Coordinates Altitude (m) Main species Coverage (%) 
15 SEP USEE IN 3678 Stipa capillata, Carex tristachya 89.00+5.13 


102°36'51.26"E 


pad Е. 3695 Stipa capillata, Anaphalis sinica 78.00+3.26 


gr ne 3702 Stipa capillata, Carex tristachya 67.00+4.15 


4% О 3721 Carex tristachya, Stipa сар Шаға 50.00+4.71 


5% ТӘ ТЫ 3725 Carex tristachya, Artemisia desertorum 38.00+3.52 
6h reer 3736 Artemisia desertorum, Agrostis matsumurae 29.00+2.20 
33°13'37.39"N 


7° 102°36/51 25"Е 3740 Artemisia desertorum, Oxytropis kansuensis 17.00+2.44 


Note: Mean+SE. 


We clipped the aboveground parts of plants by using scissors at ground level and separated 
them according to plant species. We also collected the belowground parts of plants in the 0-30 cm 
soil depth (0-10, 10-20, and 20-30 cm layers) by taking soil cores (5 cm in diameter), with three 
replicates for every layer within a quadrat. The blocks were then transported to the laboratory and 
soaked in water to remove the residual soil using a 0.5-mm sieve. The plant samples were 
oven-dried at 65°C until they reached a constant weight to determine the aboveground biomass 
(AGB) and BGB. Soil samples were also collected from the three replicate soil profiles іп the 0- 
30 cm soil depth (0-10, 10-20, and 20-30 cm layers) using 5 cm diameter soil cores and 200 cm? 
cutting rings. For the former (AGB), soil samples were handpicked to extract fine roots and other 
organic materials and were used to analyze the soil chemical properties after being air-dried and 
sieved through a 2.0-mm mesh. For the latter (ВСВ), soil samples were oven-dried at 105°C until 
they reached a constant weight, and then the soil water content (SWC) and soil bulk density 
(SBD) were measured. The soil organic carbon (SOC), soil total nitrogen (STN), soil total 
phosphorus (STP), soil available nitrogen (SAN), and soil available phosphorus (SAP) were 
determined by an external heating method (the КэСі?О? volumetric method), а vario MACRO 
cube elemental analyzer (Elementar Analysensysteme GmbH, Germany), the МаНСОз alkali 
digestion method and molybdenum antimony colorimetry, the continuous alkali-hydrolyzed 
reduction—diffusion method, and the Olsen method, respectively (Olsen et al., 1954; Bao, 2000). 


2.3 Data analysis 


The benefit of BGB in a given soil layer is expressed as the relative deviation from the mean of 
an objective observation calculated via the following equation (Bradford and Damato, 2012; Sun 
et al., 2018b): 
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pos (1) 

X max— X min 
where В; is the magnitude of the benefit for objective i; xi is the measured value of i; and Xmin and 
Xmax are the minimum and maximum values for all xi, respectively. A three-dimensional trade-off 
method was used to describe the BGB trade-offs among the three soil layers (0-10, 10-20, and 
20-30 cm) (Li et al., 2019). The vertical distance (Г) from the objective point x to the zero 
trade-off line indicates the magnitude of the benefit of BGB in each soil layer (Fig. 2). The 
trade-off was decomposed into those in the а, Б, and с directions by tracing the BGB in the soil 
layers of 0-10, 10-20, апа 20-30 cm, respectively (Fig. 2). The trade-off components in each 
direction were calculated by finding the point on the zero trade-off line using the Euclidean 
distance. 


Relative benefit of BGB, |, 


Fig. 2 Illustration of trade-offs among objective belowground biomass (ВСВ) values in the three soil layers (0- 
10, 10-20, and 20-30 cm). The zero trade-off line represents that the benefits of BGB in each layer are equal. 
BGBo-10, BGB10-20, and BGB20-30 represent the ВСВ in the 0-10, 10-20, апа 20-30 cm soil layers, respectively. 


The differences in the soil properties in different soil layers along the degradation gradients 
were calculated by using SPSS 19.0 software (SPSS Inc., Chicago, IL, USA) to determine the 
one-way analysis of variance. We used linear piecewise quantile regression analysis and linear 
regression in SigmaPlot for Windows, version 10.0 (Systat Software, Inc., Chicago, IL, USA; 
Windows version 14.0) to visualize the changes in the AGB of graminoids and forbs along the 
degradation gradients, respectively. The ggplot2 and corrgram packages in R version 3.3.2 (R 
Core Team, 2016), were employed to explore the variation in the BGB trade-off values among the 
three soil layers along the degradation gradients and to determine the correlations among the 
BGB trade-off values, soil variables, and plant community structure, respectively. 


3 Results 


3.1 Changes in soil characteristics along the degradation gradients 


SOC, STN, STP, SAN, SCOM, and SWC decreased gradually with increasing degradation 
degree, whereas no significant changes were found in SAP and SBD (Table 2). Most soil 
components such as SOC, STN, SAN, and SAP were significantly higher in the uppermost soil 
layer than in the subsoil layers. Similar patterns were found for SWC in the first and second 
degradation gradients (Table 2). 
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3.2 Variations in plant aboveground biomass (AGB) and belowground biomass (BGB) 
along the degradation gradients 


With increasing degradation degree, the AGB of graminoids decreased dramatically until the 4% 
degradation gradient and then decreased slowly, whereas the AGB of forbs changed little (Fig. 3). 
The graminoids clearly had a higher AGB than the forbs in the first degradation gradient, whereas 
the latter accounted for most of the AGB after the fourth degradation gradient. 


60 

= Graminoids (regression line) 
50 | æ Forbs (— regression line) 
404 


AGB (g/m?) 
Ww 
Ф 


1“ 2һ4 34 4th 5th 6% Th 
Degradation gradient 


Fig. 3 Changes in aboveground biomass (AGB) of graminoids and forbs along the seven degradation gradients. 
Information on the seven degradation gradients is shown in Table 1. 


The BGB clearly decreased along the degradation gradients. The value in the 0-10 cm soil 
layer (BGBo-10) decreased dramatically with increasing degradation degree (Fig. 4). By contrast, 
the ВСВ іп the 10-20 cm soil layer (BGBio-20) remained virtually unchanged along the 
degradation gradients. The ВСВ in the 20-30 cm soil layer (BGBzo0-30) tended to first decrease 
and then increase along the degradation gradients. BGB20-30 was similar їо ВСВзо о but was 
clearly lower than BGBo-10 in most of the degradation gradients, with the exception of the 7" 
gradient (Fig. 4). 
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Fig. 4 Changes of belowground biomass (ВСВ) in the three soil layers (0-10, 10-20, and 20-30 cm) along the 
seven degradation gradients 


3.3 Relationships among the BGB and soil and plant characteristics 


Figure 5 shows the variations in trade-off values of BGB іп the three soil layers (0-10, 10-20, 
and 20-30 cm) along the degradation gradients. As shown in Figure 5, the trade-off value of 
ВСВо-10 decreased remarkably with increasing degradation degree (В2-0.71, P<0.05). In the 0-10 
cm soil layer, the trade-off value of BGBo-10 exhibited significant positive correlations with SWC 
(А2=0.82), SCOM (R?=0.93), SOC (R?=0.77), STN (В?-0.76), STP (R?=0.69), and SAN 
(82=0.69), and also with the ratio of AGB of graminoids to that of forbs (і.е., AGB(G/F); 
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К2=0.66); however, it was negatively correlated with SBD (R?=0.72) at the 95% confidence level 
(Fig. 6). 
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Fig. 5 Changes in trade-off values of ВСВ in the three soil layers (0-10, 10-20, and 20-30 cm) along the seven 
degradation gradients. The gray area represents 95% confidence intervals. 
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Fig. 6 Correlations among the trade-off value of ВСВ іп the 0-10 cm soil layer (ВСВо-1), topsoil 
physical-chemical properties, and plant community structure. SBD, soil bulk density; SAP, soil available 
phosphorus; SCOM, soil compaction; SWC, soil water content; STP, soil total phosphorus; SOC, soil organic 
carbon; STN, soil total nitrogen; AGB(G/F), the ratio of AGB of graminoids to that of forbs; SAN, soil available 
nitrogen. The color depth represents correlated strength, that is, a deep color suggests a strong correlation. Blue 
and red represent positive and negative correlations, respectively. 
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4 Discussion 


41 Variations of soil properties along the degradation gradients 


In agreement with previous studies showing that degradation resulted in considerable nutrient 
loss, which then accelerated degradation (Wen et al., 2013; Pan et al., 2017; Zhang et al., 2019), 
our results indicated that soil nutrients such as SOC, STN, STP, and SAN dramatically decreased 
with increasing degradation degree (Table 2). These findings could be attributed to the following 
factors. First, the degradation-induced reductions in vegetation biomass and coverage not only 
decrease the soil organic matter (SOM) input (Wang et al., 2009), but also enhance water and soil 
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erosion, resulting in nutrient losses (Zhu et al., 2011; Dlamini et al., 2014). Second, changes in 
the microenvironment (e.g., SWC and SBD) along the degradation gradients increase the 
dissolution of SOM (Wang et al., 2008). Third, the degradation-induced turnover of species 
composition, for example, the ongoing replacement of grasses by ruderal species during grassland 
degradation, can affect the fixed soil C content (Wu et al., 2014; Han et al., 2018). Finally, 
changes in litter can have important effects on the variation of SOC (Mehta et al., 2013) and STN 
(Che et al., 2018). Specifically, degradation reduces the C and nitrogen (N) storage by changing 
the quantity and quality of litter input to the soil (Sun et al., 2018a). 

In this study, SWC decreased gradually with increasing degradation degree (Table 2), in 
agreement with numerous studies (Wei et al., 2010; Yi et al., 2012). The reason might be the 
negative effects of retrogression on the soil water-holding capacity (Fu et al., 2015) and an 
increase in evaporation owing to strong radiation due to degradation-induced plant cover 
reduction (Wang et al., 2009). Additionally, the decrease in SCOM with increasing degradation 
degree indicates that degradation causes the soil to be less cohesive by reducing the clay content 
and increasing the silt content (Yi et al., 2012). Notably, most of the soil physical-chemical 
properties were more vulnerable to degradation in the upper soil layer (0-10 cm) than in the 
subsoil layers (10-20 and 20-30 cm) (Table 1). Other works have reported similar behavior 
(Bardgett et al., 2005; Pan et al., 2017), which indicates that the impacts of degradation on the 
soil characteristics appear primarily in the upper soil layer. The reason is probably that 
degradation-induced decreases in vegetation cover and root grasp directly expose the surface soil 
to wind, rainfall, and surface flow (Pan et al., 2017). 


4.2 Effects of habitat and plant community structure оп the dynamics of BGB distribution 
patterns 


We found that the BGB distribution was depth-dependent, with more than 50% concentrated in 
the 0-10 cm soil layer, except in extremely degraded grasslands (Fig. 4). This agrees overall with 
the pattern reported by Li et al. (2011) in the alpine grassland on the Tibetan Plateau, which is 
higher than the values reported in Inner Mongolia of China (47%; Ma et al., 2008) and temperate 
zone grasslands worldwide (44%; Jackson et al., 1996). The unique environmental conditions on 
the Tibetan Plateau, such as high soil rigidity and permafrost, might contribute to the depth 
dependence (Li et al., 2011). Of most interest to our study, the trade-off value of BGB in the 0-10 
cm soil layer (BGBo-10) decreased dramatically along the degradation gradients. This finding 
could be explained by a range of direct or indirect possible mechanisms. 

The soil physical-chemical properties affect the changes in the vertical distribution of BGB. 
There is an apparent positive association between the trade-off value of BGBo-10 and SCOM in 
the topsoil (Fig. 6), which suggests that individual deeply rooted plants provide resilience against 
large gale-force winds in cohesionless soil (Vo and Johnson, 2001). Moreover, soil with higher 
SBD is generally less porous; consequently, root respiration and microbial activities, and thus root 
growth, are limited (Yan et al., 2018), which is reflected in the negative correlation between the 
trade-off value of BGBo-10 and SBDo-10 (SBD in the 0-10 cm soil layer) (Fig. 6). Additionally, 
SWC, which is tightly linked to nutrient availability, plays an important role in plant 
physiological activities (e.g., transpiration) and strongly affects plant growth (Swift et al., 2004). 
The significant positive correlation between the trade-off value of ВСВо-10 and SWCo-10 (SWC іп 
the 0-10 cm soil layer) in our study (Fig. 6) indicated that deeper-rooted plants could use water 
from the subsoil when soil water in the upper soil decreased dramatically with increasing 
degradation degree (Hu et al., 2013). 

Plants in infertile soil produce deep root systems to access nutrients from deeper soil layers 
(Wang et al., 2006; Han et al., 2011). In this study, we found significant positive correlations of 
the trade-off value of BGBo-10 with SOC, STN, STP, and SAN in the 0-10 cm soil layer (Fig. 6). 
Soil nutrients play important roles in multiple plant physiological activities: SOM provides 
energy for heterotrophic N-fixing microbes and further facilitates plant N fixation (Reed et al., 
2011; Zheng et al., 2019); soil N has important effects on the leaf N and plant photosynthesis 
(Yue et al., 2019); low soil P can greatly limit plant growth and functioning (Quesada et al., 2009; 
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Sun et al., 2019b). Consequently, soil nutrient changes can result in species turnover by affecting 
plant facilitation and competitive exclusion (Marriott et al., 2002; Guo, 2003). 

Dynamics in the plant community structure directly affect the BGB vertical distribution 
patterns. Our results showed a dramatic decrease in AGB of graminoids and a slight increase in 
AGB of forbs with increasing degradation degree (Fig. 3). This finding suggests that plant 
communities tend to shift from graminoids to forbs after degradation occurs, showing successive 
vegetation reversal with the degradation of the plant community structure (Wang et al., 2009; Wu 
et al., 2009). This result is mainly because that the main cause of grassland degradation in our 
study site is free grazing, and selective foraging decreases the proportion of palatable graminoids. 
Moreover, forbs show a stronger competitive capacity than graminoids in the increasingly hostile 
environment resulting from increasing degradation degree (Sun et al., 2009; Han et al., 2011). 
Graminoids have shallow roots that extend horizontally in the surface soil, whereas forbs produce 
deep axial roots (Wang et al., 2004), which contribute directly to the reduction in the root biomass 
in the topsoil, as indicated by the significant positive correlation between the trade-off value of 
BGBo-10 and AGB(G/F) (Fig. 6). 


5 Conclusions 


In the current study, an alpine meadow with sequential degrees of degradation at the Zoige 
Plateau was selected as a typical degraded grassland to analyze the variations in plant BGB 
distribution under degradation. The trade-off value of BGB in the soil layer of 0-10 cm clearly 
decreased with increasing degradation degree and exhibited significant positive correlations with 
SWC, SCOM, SOC, STN, and STP. Our findings suggest that alpine grassland degradation 
dramatically reduces belowground net primary production in the topsoil owing to 
degradation-induced deterioration of soil conditions and changes in the plant community 
structure. These results can provide valuable information for scientific assessment of the effects 
of degradation on the soil C pool distribution in alpine meadows on the Tibetan Plateau and 
similar areas worldwide. Larger-scale transect surveys should be used in future studies to further 
explore the general response of root biomass distribution to grassland degradation. 
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